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Abstract 

Finite Element simulations of rubbers and biological soft tissue 
usually assume that the material being deformed is slightly compress- 
ible. It is shown here that in shearing deformations the corresponding 
normal stress distribution can exhibit extreme sensitivity to changes in 
Poisson's ratio. These changes can even lead to a reversal of the usual 
Poynting effect. Therefore the usual practice of arbitrarily choosing a 
value of Poisson's ratio when numerically modelling rubbers and soft 
tissue will, almost certainly, lead to a significant difference between 
the simulated and actual normal stresses in a sheared block because 
of the difference between the assumed and actual value of Poisson's 
ratio. The worrying conclusion is that simulations based on arbitrar- 
ily specifying Poisson's ratio close to 1/2 cannot accurately predict 
the normal stress distribution even for the simplest of shearing defor- 
mations. It is shown analytically that this sensitivity is due to the 
small volume changes which inevitably accompany all deformations of 
rubber-like materials. To minimise these effects, great care should be 
exercised to accurately determine Poisson's ratio before simulations 
begin. 
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1 Introduction 



When solid rubber-like materials are assumed to be incompressible so that 
only isochoric (i.e. volume preserving) deformations are allowed, then the 
theory of non-linear elasticity is known to be extremely effective in providing 
analytical predictions for their mechanical response when both the original 
geometry and boundary conditions are simple. These analytical solutions 
include shear, torsion, bending, straightening and inflation (see, e.g., Rivlin 
[2], Green and Zerna |3], or Ogden |4]). To analyze more complex prob- 
lems, use of Finite Element (FE) simulations is usually necessary. However, 
in most commercial finite element codes, incompressibility is not assumed 
ab initio in order to prevent element locking (see, e.g., |5l El [^). Thus in 
FE analyses, slight compressibility is usually assumed and the solution for 
fully incompressible materials is obtained in a limit process if necessary (see, 
e.g., [H El El Hn])- Indeed, because all solid rubbers are to some extent 
compressible, models for slight compressibility have been investigated inde- 
pendently of any FE considerations (see, e.g., Ogden [1], Horgan and Murphy 
[El [HI Uni US] and references cited therein) . 

Rubber components are often subjected to shearing deformations in ap- 
plications and biological soft tissue is often sheared in vivo (see, for example, 
LeGrice et al. [T7| for a discussion on shear strain in the left ventricular my- 
ocardium or Horgan and Murphy [18] on the simple shearing of tissues). It is 
therefore essential that accurate, reliable and efficient models be available to 
predict shear behaviour. In a recent paper, Gent et al. [1] conducted some 
numerical simulations using the commercial code ABAQUS of the simplest 
possible shearing deformation, where one face of a block, modelled as a com- 
pressible neo-Hookean solid, is displaced relative to the parallel face. This 
shearing deformation will be called here experimental simple shear. One set 
of their results is puzzling: in their Figures 12 and 13, an infinitesimal change 
in Poisson's ratio (z/) leads to a large change in the predicted values of the 
normal stress components, with, in some cases, a small change in Poisson's 
ratio leading to a change in the sign in the components, leading effectively to 
an 'inverse' Poynting effect. These numerical results have been confirmed by 
our own numerical simulations (§2) for which the adopted procedure mirrors, 
as closely as possible, the earlier work. 

In §3 this extreme sensitivity of the normal stresses to changes in v is 
explained by combining a realistic mathematical model of experimental sim- 
ple shear that incorporates the possibility of small volume changes with the 
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constitutive models of slight compressibility used in FE simulations. It will 
be shown that the predicted normal stresses for the slightly compressible 
neo-Hookean material obtained from this combination display exactly the 
sensitivity to changes in the Poisson ratio as was first observed by Gent et 
al. pp. The shearing strains considered here are only moderate in range, 
with a maximum shear strain of 100% being imposed. It is likely that the 
extreme sensitivity noted here is more pronounced for larger strains. 

We suggest that the results discussed above have significant implications 
for the numerical simulation of materials traditionally considered as being 
incompressible. We note that in practice, when simulating these materials, 
experimental values of the corresponding Poisson ratio are not known. In- 
stead, most engineers adopt the pragmatic solution of specifying instead a 
value of Poisson's ratio 'close' to 1/2, assuming implicitly that the specific 
value of Poisson's ratio adopted will not have a significant impact on the nu- 
merical results and, in particular, that small percentage changes in the value 
of Poisson's ratio result in corresponding small percentage changes in the 
stress distribution. Most simulations of biological soft tissue, for example, 
are conducted on this basis. However, the simulations of Gent et al. [1] and 
our own numerical experiments (§2) all suggest that, with the models cur- 
rently used, a small change in the value of Poisson's ratio chosen (say from 
0.499 to 0.495) will result in significant changes in the stress distribution. 
There is almost certainly a difference between the assumed and actual values 
of Poisson's ratio and this should generally result in a significant difference 
between predicted and actual stress distributions. Our conclusion is that an 
accurate determination of Poisson's ratio is essential if FE simulations are 
to yield accurate normal stresses in shearing deformations (See also Horgan 
and Murphy [T^] for a discussion of the results of Gent et al. [I] relevant to 
the present paper.) 

It is finally noted that the seemingly obvious solution to this extreme 
sensitivity to Poisson's ratio of simply simulating perfect incompressibility is 
not a valid approach. No material is perfectly incompressible and assuming 
this idealisation (denoted in our results by z/ = 1/2) will result in a difference 
between the idealised value assumed and the actual value of Poisson's ratio. 
This difference will lead to a significant percentage error in the predicted 
normal stresses. 
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2 Modelling slight compressibility 



Many different constitutive models have been proposed to reffect slight devi- 
ations from incompressibility on assuming that the material is homogeneous, 
isotropic and hyperelastic. We refer to the recent papers [131 EH [ISl ES] for 
background, references to the pertinent literature and, very importantly, a 
summary of the fit of these models, and, in particular, the models used here, 
with experimental data. Here we briefly describe the usual form of such 
slightly compressible (or almost incompressible) strain-energy functions used 
in the FE simulations. 

Assume that an incompressible strain-energy function W = (Ai, A2, A3), 
where ^ is a symmetric function of the principal stretches Aj, has been ob- 
tained that gives an excellent fit to the experimental data collected from 
some set of material characterization experiments. Then the corresponding 
slightly compressible form implemented in FE models is usually of the form 

W^ = ^(Ai,A2,A3)+F(23), (1) 

where 

is = A1A2A3, Ai = Ai/4^^, (2) 

F is some specified polynomial function of its indicated argument, and the 
A's are the so-called 'deviatoric stretches' (see the ABAQUS [H] and ANSYS 
[12] manuals). 

On assuming that (1, 1, 1) = to ensure zero strain-energy in the ref- 
erence configuration for the incompressible material, it follows from ([T]) that 
we need to assume that 

F(l) = 0, (3) 

for the same property to hold for the slightly compressible material. To 
ensure zero stress in the reference configuration, it is sufficient to require 
that 

F'(l) = 0. (4) 

Finally, to ensure compatibility with the classical linear form of the elas- 
tic strain-energy function on restriction to infinitesimal deformations, it is 
required that 
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where n and k are the infinitesimal shear and bulk moduli, respectively. 
Recall that they are connected to Poisson's ratio u through 

= 2(1 + z/)/[3(l-2z/)]. (6) 

Motivated by these relations, it is now assumed that the ijj term in ([T]) 
is proportional to /i and that F is proportional to k. On setting is = 1 in 
([T]), one recovers the original strain-energy function for the corresponding 
incompressible material. It is important to emphasise that the main moti- 
vation for strain-energy functions of the form ([T]) seems to be mathematical 
convenience: a simple additive term in is included to account for the com- 
pressibility so that, on restriction to infinitesimal deformations, and F are 
proportional to the shear modulus fi and the bulk modulus k, respectively 
(see e.g., Horgan and Murphy [16] for further details). 

Many FE codes are based on models of the form ([T]). The initial moduli 
are usually assumed to be such that k/^ 3> 1 in order to model slight in- 
compressibility. In the variational formulation, the bulk modulus k, is called 
the penalty number (see, e.g., [5l |20l [7]) and is chosen to be large. From 
IQ, we find that for u = 0.49, 0.499, 0.4999, we have n/fx ~ 50, 500, 5000, re- 
spectively. The default setting in Abaqus/Explicit is quite low, as it assumes 
that = 20, corresponding to a Poisson ratio of z/ ~ 0.475. 

Equivalently, the model can be written in terms of the Cauchy-Green 
strain invariants, defined in terms of the principal stretches Aj, i = 1,2, 3, as 
follows, 

h = \l + \l + A3, I2 = XlXl + A^Ag + A2A3, J3 = = A^AgAg. (7) 

The corresponding strain-energy function that is used in FE simulations to 
model almost incompressible behavior then has the equivalent form 

W = ^ {h,h) + F{iy^), where h = h/ll'\ h = h/lT- (8) 

Further details can be found, for example, in Crisfield [TO] . 

The general representation of the stress-strain relations for compressible, 
non-linear elasticity can be found, for example, in Atkin and Fox [21] and is 
given by 

T = «oi" + aiB + aa-B^, (9) 

where 

^.mdW 2 fdW ^ dW\ 2 dW 
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Here T is the Cauchy stress, / is the identity tensor, and B is the left Cauchy- 
Green strain tensor, defined in terms of the deformation gradient tensor F 
through the relation B = FF'^ . The stress-strain relation for the almost 
incompressible material, ([s]), therefore has the form (|9]) with the following 
coefficients: 

2 J dip 4 jdip 
_ 2 (dip -dip 

The stress-strain relation for almost incompressible materials can thus be 
put in the form 

T = F'{i^)I + aQl + aiB + a2B'^, (12) 



with ai, a2 given in (11)2,3 and do = ao — F'{is), i.e. 

2 dip A dip 

Thus the stress is always the sum of two parts: a hydrostatic term propor- 
tional to the volumetric function -^(^3), and therefore of the order of the 
infinitesimal bulk modulus n, and another term of the order of the infinites- 
imal shear modulus /x. 

The numerical simulations of Gent et al. [1] are based on the slightly 
compressible neo-Hookean specification of this strain-energy function, readily 
implemented in ABAQUS [TT]. It has the form 

W = ^{h-?,)+F (^3) = ^ (A? + + - 3) + F (.3) , (14) 

where the constitutive volumetric function F is the quadratic form 

n^-^) = ^ (^3 - 1)' , (15) 

which is the simplest form to satisfy (|3]), (|4]), and For this material the 
constitutive coefficients simplify to 

"o = -r7H7^' "1 = 757^' "2 = 0, (16) 

<J-'3 -'3 
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v=Q.47 



Figure 1: Normal stress variations along the horizontal center line of a 
sheared block of a slightly compressible neo-Hookean material (amount of 
shear: 1, block's aspect ratio: 10). The Cauchy stress components Tn (left) 
and T22 (right) are highly sensitive to small changes in the values of Poisson's 
ratio p = 0.5, 0.499, 0.495, 0.49, 0.48, 0.47. 
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Figure 2: Finite element simulation of the simple shear of a block with the 
same dimensions and characteristics as the block of Gent et al. [T] (amount 
of shear: 1, block's aspect ratio: 10). Here Poisson's ratio is z/ = 0.5 (perfect 
incompressibility) and the normal stress T22 is tensile throughout the block, 
except in small areas near the compressed corners. 



and the following stress-strain relation is obtained: 

T = ^{li^' - 1)/ + -g^ - Ihl) . (17) 
-'3 

Predictions based on the neo-Hookean model have been found to yield a 
good fit with data from a variety of experiments on rubber-like materials for 
small to moderate strains. Since moderate strains are only considered here, 
the results obtained in our analyses are therefore likely to be obtained from 
simulations based on other forms of the deviatoric part of the strain-energy 
function. 



The form ( 14 ) was assumed in our repeating, and furthering, of the numer- 
ical simulations of Gent et al. [Ij. We entered their physical and geometrical 
parameters into ABAQUS 6.7 to try to replicate their results for a sheared 
neo-Hoohean block. We obtained the stresses and strains throughout the 
thickness. We varied the mesh resolution, the time step duration and the 
number of increments used and found virtually no difference with their sim- 
ulations. Figure [2] shows an example of the stress field distribution in the 
sheared block, displaying clearly an almost homogeneous field except near 
the edges. Further pictures obtained for the stress and strain distributions 
are identical to those of Gent et al. P, see Figure [l] 
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We also carried out two re-meshes in addition to the original mesh pQ 
of 20 (thickness) by 100 (length), i.e. 2,000 rectangular elements. The first 
re-mesh was a square mesh of 20 x 200 (i.e. 4,000 elements) and the second 
re-mesh was a square mesh of 40 x 400 (i.e. 16,000 elements). Each new 
mesh was run for the two following values of Poisson's ratio: v = 0.49 and 
u = 0.5. These re-meshes had no effect on the original results of Gent et al. 
[1]. Similarly, we varied the time step and incrementation, from their original 
[1] values of time step = 4 seconds, 400 increments (giving a shearing rate of 
2 mm/s). Thus we halved (time=2, increments=200, shearing rate=4mm/s) 
and then doubled (time=8, increments=800, shearing rate=lmm/s) the time 
steps and increments. These time changes also has no effect. 

In conclusion, we found, just like Gent et al. [1] , that for the compressible 
neo-Hookean model (14)-(15) infinitesimal changes in the value of Poisson's 
ratio can lead to finite changes in the simulated stress distribution. Our 
simulations suggest that this problem is especially acute for z/ G [0.495,0.5] 
and that that the problem of extreme sensitivity of the normal stresses to 
Poisson's ratio increases the closer the value of u is to 1/2. These features 
are evident in Figure [T] and, perhaps more emphatically, in a comparison 
of Figures [2] and [3j In Figure |2| where perfect incompressibility has been 
assumed, the stress is almost homogeneous, except for more complicated be- 
haviour near the corners, and tensile throughout a block with aspect ratio 
1:10 sheared by a unity amount. Figure [3] displays the normal stress compo- 
nent T22 distribution with Poisson's ratio u = 0.495. In contrast with Figure 
|2} the stress is no longer homogenous: it can be seen there that T22 is ten- 
sile and non-uniform near the slanted faces, but is compressive towards the 
centre of the specimen. 

It is natural to see if the problems in simulating slightly compressible 
materials are reflected in analytical approaches to simple shear because their 
specific cause might then be isolated. To investigate these problems analyti- 
cally, the classical mathematical approximation to experimental simple shear 
is generalised in the next section and the corresponding stress distribution 
obtained for both the general slightly compressible model (12) and for the 
slightly compressible neo-Hookean model (17) used in FE simulations. 
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Figure 3: Finite element simulation of the simple shear of a block with the 
same dimensions and characteristics as the block of Gent et al. [Ij (amount 
of shear: 1, block's aspect ratio: 10). Here Poisson's ratio is z/ = 0.495 (slight 
compressibility) and the normal stress T22 is uniform and compressive in a 
large part of the center of the block. 



3 Analytical approaches to simple shear 

The classical mathematical model of simple shear assumes that a block of 
rubber is deformed according to 

a;i = Xi+7X2, X2 = X2, X3 = X3, (18) 

where (Xi, X2, X3) and (xi, X2, X3) denote the Cartesian coordinates of a typ- 
ical particle before and after deformation, respectively, and 7 is an arbitrary 
dimensionless constant called the amount of shear. Within the context of 
nonlinear elasticity, this problem was first considered by Rivlin [2j and has 
been much discussed and analyzed since. The reference works of Green and 
Zerna [3], Truesdell and Noll [22], and Ogden P], for example, all consider 
various aspects of this problem. 

Intuitively it seems reasonable that this homogeneous deformation should 
accurately model experimental simple shear, especially within the bulk of the 
material and for modest amounts of shear. The corresponding in-plane stress 
distribution for the general strain-energy function adopted in FE simulations 



to model almost incompressible rubbers is easily obtained from (11), (12) 
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and is given by 

Til = ao + ai (l + 7^) + as (l + + 7^) , 
T22 = 60 + «! + "2 (1 + 7^) , 

T12 = ai7 + ^27 (2 + 7^) . (19) 



Modelling experimental simple shear using the intuitive deformation ( 18 ) 
therefore results in the in-plane stress field for almost incompressible mate- 
rials of the form ([T]) being independent of the volumetric function F. This 
is not surprising given that classical simple shear is an isochoric deformation 
(zs = 1). Therefore, for a fixed value of /x, the in-plane stress field is inde- 
pendent of the bulk modulus /t, or equivalently, of the Poisson ratio v. This 
is in conflict with the numerical predictions discussed earlier. 

This conflict can be explained by the fact that using the isochoric ap- 



proximation (18) to model experimental simple shear is too restrictive and 
does not capture an essential feature of every deformation of rubbers, i.e., 
that an infinitesimal volume change always occurs and therefore occurs in 
experimental simple shear in particular. To get an idea of the order of the 
change, we rely on the classical experimental data of Penn • In a series of 
elegant experiments on peroxide vulcanizates of natural gum rubber using a 
dilatometer technique, Penn measured the volume change in simple tension 
for axial stretches of the order of 2. The volume change, calculated as — 1, 
was found to be of the order of 0.0001. No comparable data exist for simple 
shear. We now investigate the effect, if any, of a change in volume. 

We model the infinitesimal volume change with the following simple gen- 



eralization of the classical approach (18) 



xi = (l+ei)Xi+7(l+e2)X2, X2 = (l + e2)X2, X3 = (l + e3)X3, 

(20) 

where the e's are infinitesimal quantities. This deformation seems intuitively 
to capture the main features of the deformation of a block when one face 
is displaced relative to a parallel face. It corresponds to a simple shear of 
amount 7 superposed on a triaxial extension, where the e's are the infinites- 
imal elongations. 

Retaining only the first-order terms in the e's gives Z3 ~ 1 + e, where 
e = ei + e2 + e3 is the volume change. It follows that, to the first-order in 



e, the constitutive function F'{i^) occurring in the stress-strain relation (12) 
for almost incompressible materials can be approximated as 

F\iz) ^ ^^e, (21) 
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where we used Q and 

The corresponding in-plane stresses are given by Rajagopal and Wineman 

as 

Tn = Ke + ao + tti(A? + + aa^A? + -f^Xj^ + -f^X% 
T22 = Ke + ao + aiXl + 0:2(1 + 7^)A2, 

T12 = ai7A^ + «27A^[A? + (1 + 7')A^], (22) 

where Aj = 1 + are the stretch ratios (here it is understood that their 
powers are hnearized with respect to the e's). Compare these general ex- 



pressions with those for simple shear given in (19): the shear stress T12 is 
again independent of the bulk modulus and volume change, but now both 
normal stresses Tn and T22 contain a new ne term. This is in agreement with 
the effects observed in the numerical simulations, showing that the normal 
stresses, not the shear stress, are sensitive to volume changes (see Figures [l] 
andga)). 

Now we try to mirror the numerical simulations: we specialize the stress- 



strain relation to that of the slightly compressible neo-Hookean solid (17); we 
impose a parallel displacement for the moving face with respect to the fixed 
face (so that A2 = 1); and we impose a two-dimensional formulation (so that 
A3 = 1). From [21] we find that the first principal invariant is Ii = 3 + 7^ + 2e, 
where we linearized with respect to e = ei. For this specialisation the normal 
stress component T22, for example, is easily computed and we find that 

T22 7^ / 2z/ 57 



where we used the connection ^ between the bulk modulus, the shear mod- 
ulus and Poisson's ratio. To compare this stress with the simulated stress 
given in Figure |2| the volume change e = 23 — 1 is needed. In Figure |4]^b) is, 
as measured by ABAQUS, is plotted as a function of the amount of shear 7 
for different values of Poisson's ratio: u = 0.499, 0.495, 0.49. In Table 1, we 
collected in particular results obtained when the amount of shear is 7 = 1.0, 
as in the simulations. The volume change ^3 — 1 is displayed in the second 



column and the normalized Cauchy stress component T22//i, found from (23 ) 
in the third column. 



It follows that (23) predicts normal stresses that are very close to the 
average values of T22 predicted by the numerical simulations given in Figure 
^h). Therefore infinitesimal changes in Poisson's ratio close to the limiting 
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Figure 4: On the left: Shear stress variation along the horizontal center line 
of a sheared block of a slightly compressible neo-Hookean material (amount 
of shear: 1, block's aspect ratio: 10). The Cauchy stress component T12 
is not very sensitive to small changes in the values of Poisson's ratio u = 
0.5, 0.499, 0.495, 0.49. On the right: Volume change during the large shear of 
a slightly compressible block of a neo-Hookean material with an aspect ratio 
of 10. 
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V Z3 - 1 



0.499 0.0009 0.116 
0.495 0.0029 -0.045 
0.490 0.0044 -0.115 

Table 1: Variation of volume change with Poisson's ratio 



value of V = 1/2 may indeed result in dramatic changes in the value for the 
normal stress, even leading to an inverse Poynting effect, due to the volume 
changes which accompany every deformation. 



4 Conclusion 

It has been shown that the usual practice of arbitrarily choosing a value of 
Poisson's ratio when numerically modeling rubbers results in widely differ- 
ent normal stress distributions when the third decimal of z/ G [0.495, 0.5] is 
modified. The sign of the Poynting effect may even be reversed. It is thus 
essential for the accurate modeling of rubbers to experimentally determine 
Poisson's ratio to the utmost precision. This is not as hard as it seems, es- 
pecially when the measure of v relies on ultrasonic wave speeds evaluation. 
There, the contrast between the speed of a compression wave vl and the 
speed of a shear wave vt is so large for nearly incompressible solids, that the 
connection 

^=2—2 f' 24 

gives excellent precision. For instance. Wood and Martin find u = 
0.49986 using acoustic waves. Similarly, Gennisson et al. |26] measure the 
following body wave speeds for an Agar- Agar gel: vl = 1500 m/s and f t = 1 



m/s, giving u = 0.49999978, according to (24). 

The issues raised here ought to be addressed by anyone wishing to use 
a commercial FE code for large deformation behavior. For instance the 
Abaqus/Explicit manual [TT] indicates that the default setting for shght com- 
pressible models is u = 0.475, a value which allows volume changes sufficient 
to inverse the Poynting effect according to Figure 1. It goes on to state 
that if we "are defining the compressibility rather than accepting the default 
value, an upper limit of 100 is suggested for the ratio of which corre- 
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sponds to 1/ ~ 0.495, also within the range of extreme sensitivity in shear. 

Furthermore, the manual goes on to say that slight compressibility "does not 
warrant special attention for plane stress" , a statement clearly at odds with 
our findings. 
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